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Description of ^-distributions at RHIC energies 
in terms of a stochastic model 
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To explain jy-distributions at RHIC energies we consider the Ornstein-Uhlenbeck process. 
To account for hadrons produced in the central region, we assume existence of third source 
located there (y « 0) in addition to two sources located at the beam and target rapidities 
(±j/ m ax = ±ln [y/SNN /ttin])- This results in better x 2 /n.d.f. than those for only two sources 
when analysing data. 

We are interested in analyses of ^-distributions^ ^ at RHIC energies by means 
of stochastic process, in particular the Ornstein-Uhlenbeck (O-U) process.^® Our 
previous formulation was a simplified approach, because we have sought an empirical 
formula. Through those studies, we have found that 1) ^-distributions should be a 
sum of two or more Gaussian distributions, and 2) a z r = r//ry rms scaling holds among 
various centrality cuts. We have used the following Fokker-Planck equation 
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where y, 7 and a 2 are the rapidity, the frictional coefficient and the variance, respec- 
tively. The solution with the initial condition P(y, 0) = 0.5[<5(y + y max ) + <5(y — y ma x)] 
is given as 
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where V?(t) = {a 2 /2 7 )p with p = (1 - e~ 2 ^). See Refs. 5) and©. 

In this paper we investigate a more realistic formulation, namely we take into 
account the contribution from the central region at y OP For our aim, first of all, 
we adopt the two-step processes by O-U process shown in Fig. ^ 

For the third source we assume the following y distribution 
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where Po(y, 0) = 5(y) and U 2 (i) = (ao/2 r y)p. According to the physical picture 
mentioned above we have therefore following formula for the normalized distribution 
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(J dn/dn ■ dn = 1) in the //-rapidity space 
dn 



dn 



J(rj, m/p t ) x 77 x Eq. Q + J x (1 - 77) x Eq. © 
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where the Jacobian factor J(r), m/pt) = cosh77/y 1 + (m/pt) 2 + sinh 2 77 and Jo = 
J(t7, mo/pt). rt is the ratio of the weight factors among the two sources (ci) at 
iymax and the third one (co); 77 = 2ci/(2ci + Co). 

In concrete analyses, the masses of hadrons in the beam and target nuclei 
should be larger than those of the central region, because of the richness of baryons; 
(m/pt) > (mo/pt), (m/pt) > 5(m/p t ), (mo/pt) > S{mo/pt)- The evolution parame- 
ter p is determined at the minimum values of x 2 ' s > as h 1 t ne previous treatments 

Our analyses of data by PHOBOS Collaboration^ at ^snn = 130 GeV and 
200 GeV are shown in Fig. [21 and Table [IJ As one can see there, a reasonable set of 
weights, which satisfy our criteria, is c\/2 : cq : ci/2 = 1 : 6 : 1 ~ 1 : 10 : 1. 

Moreover, to confirm the z r = rj/r] rrns (7/ rms = y {rj 2 )) scaling we compare our 
theoretical formula 
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with data including 5 centrality cuts. Explanations of the z r scaling by means of 
Eq. © shown in Fig. 01 seem to be excellent. 

Hereafter, we analyzed data by BRAHMS Collaboration 2 ' at 200 GeV using 
Eqs. @ and ©. In this case our criterion mentioned above cannot be applied, 
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Table I. Weights (ci/2 : Co : ci/2 = 1 : Cq : 1) dependence of x 2 's values. 





CO 


2 


4 


6 


7 


8 


10 


PHOBOS 


mo/pt 


0.29±0.91 


0.33±0.56 


0.46±0.11 


0.51±0.28 


0.55±0.25 


0.63±0.20 


130 GeV 


m/pt 


0.55±0.18 


0.54±0.20 


0.68±0.15 


0.72±0.20 


0.76±0.21 


0.83±0.22 


0-6 % 


X 2 /n.d.f. 


0.79/47 


0.71/47 


0.64/47 


0.62/47 


0.62/47 


0.65/47 


PHOBOS 


m/pt 


0.22±1.02 


0.36±0.47 


0.47±0.34 


0.51±0.29 


0.55±0.27 


0.61±0.23 


200 GeV 


m/pt 


0.45±0.16 


0.54±0.24 


0.69±0.28 


0.58±0.26 


0.62±0.24 


0.69±0.33 


0-6 % 


X 2 /n.d.f. 


0.63/47 


0.68/47 


0.76/47 


0.81/47 


0.88/47 


1.02/47 




Fig. 3. z r scaling by Eq. Q. (a) xV n - d - f - = 21.2/318. (b) xV n - d - f - = H-5/318. Dotted lines 
mean the error bars. 



as seen in Table ITTl because of fluctuation in dn/drj and restricted data points in 
fragmentation regions. Therefore we have assumed the evolution parameter p = 0.75 
taken from Fig. Efb) which results in reasonable value of x 2 ' s a t ci/2 : Co : c\j1 = 
1 : 20 : 1. The z r scaling shown in Fig. 0] holds with ci/2 : cq : c\/2 = 1 : 10 : 1 in a 
sense of the averaged centrality cuts. 

In conclusion, the two-step processes by O-U process is available for the descrip- 
tion of ^-distributions at RHIC energies. 
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Table II. Same as Table [Q but data by BRAHMS Collaboration at 200 GeV. 



/ 

Co 


5 


p = free 
10 


20 


p = 0.75 (fixed) 
5 10 20 


BRAHMS 
200 GeV 
0-5 % 


mo/pt 
m/pt 
P 

X 2 /n.d.f. 


0.57±0.33 
6.2±2.9 

1.00±0.80 
4.70/30 


0.51±0.18 
7.1±2.7 

1.00±0.84 
4.56/30 


0.56±0.16 
16.7±25.0 
1.00±0.03 
3.89/30 


0.29±0.38 0.37±0.46 0.55±0.24 
0.00±2.15 0.35±1.57 0.78±0.66 
0.75 0.75 0.75 
4.94/31 4.70/31 4.37/31 
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Fig. 4. (a) ci/2 : c : ci/2 = 1 : 20 : 1, X 2 = 4.68/30, y max = 5.36 (fixed), p = 0.75 (fixed), 
Vo = 6.42 ± 0.95, Vi = 0.79 ± 0.23. (b) z r scaling, ci/2 : Co : ci/2 = 1 : 10 : 1, p = 0.95 ± 0.06, 
X 2 /n.d.f. = 29.8/186. Dotted lines mean the error bars. 
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